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Abstract 

We present a unified description of the relativistic irNN and ^yirNN sys- 
tems where the strong interactions are described non-perturbatively by four- 
dimensional integral equations. A feature of our approach is that the photon is 
coupled in all possible ways to the strong interaction contributions. Thus the 
hadronic processes NN—*NN, NN -^nd, Trd^ird, etc. and corresponding elec- 
tromagnetic processes NN — > jNN, jd — > NN, jd — > nd, ed — ► ed, ed — ► eVti, 
etc., are described simultaneously within the one model of strong interactions. 
Our formulation obeys two and three-body unitarity, and as photons are coupled 
everywhere in the strong interaction model, gauge invariance is implemented 
in the way prescribed by quantum field theory. Our formulation is also free 
from the overcounting and undercounting problems plaguing four-dimensional 
descriptions of irNN-hke systems. The unified description is achieved through 
the use of the recently introduced gauging of equations method. 

I. INTRODUCTION 

Recently we have introduced a method for incorporating an external electromagnetic field 
into any model of hadrons whose strong interactions are described through the solution of 
integral equations . The method involves the gauging of the integral equations themselves, 
and results in electromagnetic amplitudes where an external photon is coupled to every part of 
every strong interaction diagram in the model. Current conservation is therefore implemented in 
the theoretically correct fashion, i.e. as prescribed by quantum field theory. Initially we applied 
our gauging procedure to the relativistic three-nucleon problem whose strong interactions are 
described by standard four-dimensional three-body integral equations fUJl]. More recently we 
used the same method to gauge the three-dimensional spectator equation for a system of three- 
nucleons H|J. 

Here we apply our gauging procedure to the more complicated case of the relativistic tiNN 
system whose four-dimensional integral equations have only recently been derived |||7| .[] These 
equations obey both 2-body (NN) and three-body (ttNN) unitarity and simultaneously describe 
all the strong interaction processes of the nNN system, including the reactions NN —>■ NN, 
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FIG. 1. Example of overcounting in NN^nd. (a) The NN^nd Feynman diagram where dark 
circles represent all possible contributions, (b) One of the contributions included in (a), (c) Another 
way of drawing diagram (b) showing that this term corresponds to an overdressing of the deuteron 
vertex. 

NN — > nd, NN — > nNN, nd — > nd and nd — > nNN. After gauging these nNN equations 
we obtain gauge invariant expressions for all possible electromagnetic processes of the nNN 
system, e.g. pion photoproduction 7c/ — > nd and 7c? — > nNN, pion electroproduction ed — > e'nd 
and ed — > e'nNN, and because pion absorption is taken into account explicitly, we also obtain 
gauge invariant expressions for processes like deuteron photodisintegration 7c? — > NN and 
Bremsstrahlung NN — > 7AA that are valid even at energies above pion production threshold. 
Included in the electromagnetic processes described by our model is the especially interesting 
case of elastic electron-deuteron scattering. Here the deuteron is described within our model 
as a bound state of the nNN system; thus, after gauging, our model provides a rich description 
of the electromagnetic form factors of the deuteron with all possible meson exchange currents 
taken into account. The four-dimensional integral equations for the nNN system are a good deal 
more complicated than the corresponding equations describing a strictly three-body system like 
that of three nucleons. Firstly the nNN equations take into account pion absorption so there is 
coupling to the two-body AW channel. A consequence of this in the four-dimensional approach 
is that certain three-body forces must be retained both to avoid the undercounting of diagrams 
and to ensure that the equations satisfy two- and three-body unitarity ||. Also, unlike in the 
NNN case, not all three-body reducible diagrams of the nNN system have unique three-body 
cuts. This leads to an overcounting of diagrams problem in the formulation of four-dimensional 
nNN equations that has only recently been solved ||. One example of the type of overcounting 
encountered is given in Fig. |l|. The newly derived nNN equations have the feature that their 
kernels contain explicit subtraction terms which eliminate all such overcounting. 

The presence of pion absorption and the overcounting problem in the nNN system would 
normally make the task of formulating a description of the jnNN system an especially difficult 
one. To put this difficulty into perspective, we recall our experience with gauging the four- 
dimensional three- nucleon system |T|,0j. In this case the strong interaction integral equations 
have exactly the same form as the three-body equations of quantum mechanics (QM) (except 
of course that they are four-dimensional rather than three-dimensional). Then after applying 
our gauging procedure we obtained an electromagnetic current that is again just of the same 
form as the three-nucleon current in QM; namely, a simple sum of one and two-body currents. 
However, in contrast to the QM case, the two-body current obtained in the four-dimensional 
formulation is modified by a subtraction term whose presence is necessary to stop the over- 
counting of diagrams. Thus the gauging of equations method led us to a simple prescription 
for the electromagnetic current of three nucleons. The situation for nNN is quite different as 
the strong interaction equations for this system cannot be cast into a QM form. Thus there 
is no corresponding simple prescription for the nNN electromagnetic current and the gauging 
procedure itself becomes effectively the only way to specify this current. On the other hand, our 
gauging procedure is extremely simple, and by gauging the nNN equations with subtraction 
terms included, one easily constructs equations for the jnNN system without encountering any 
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further difficulties and with all overcounting problems being taken care of automatically. In 
this way we obtain a unified description of the nNN and jnNN systems. 



II. FOUR-DIMENSIONAL tvNN EQUATIONS 

The first attempts to formulate few-body equations using relativistic quantum field theory 
were made already in the early 1960's |3-nJ. Both such general formulations and ones more 



specific to the irNN system have been pursued until the present time ITTHTH]. Yet all these 



attempts have had theoretical inconsistencies, including the undercounting and overcounting 
problems discussed above. We have recently overcome these problems and derived new consis- 
tent nNN equations by using a method where in cases like that of Figs. |l|(b) and (c) where the 
rightmost ttNN cut is not unique, one of the rightmost irNN vertices is "pulled out" further 
to the right, in this way defining a unique ttNN cut ||. The same ttNN equations were later 
derived in Ref. J7j where a method based on Taylor's classification of diagrams |8|Jl4}| was used. 
In this section we would simply like to restate these equations but in a form that is particularly 
convenient for gauging. In this section we would simply like to restate these equations but in 
a form that is particularly convenient for gauging. 



A. Distinguishable nucleon case 

Initially we would like to consider the case where the two nucleons are treated as distin- 
guishable particles. Not only does this avoid taking into account the symmetry factors and 
related complications due to the identity of the two nucleons, but it is also of practical interest 
in itself, as for example in modelling the imp system without using isospin symmetry. 

We follow the usual convention and refer to the two nucleons as particles 1 and 2, and the 
pion as particle 3. We also use A = 1 or 2 to label the channel where nucleon A and the pion 
form a two-particle subsystem with the other nucleon being a spectator, and A = 3 to label the 
channel where the two nucleons form a subsystem with the pion being a spectator. Although 
we are concerned only with those physical processes having at most one pion in either initial 
or final state, we do allow multiple-pion intermediate states. These intermediate state pions 
may be taken to be either distinguishable or indistinguishable without affecting the formulation 
below in any essential way. 

Although the derivation of the irNN equations given in Ref. || neglected all connected 
diagrams that are both NN- and it AW-irreducible in the processes NN — ► NN and NN <-» irNN, 
in our approach these are easily included and do not complicate the original irNN equations 
in any essential way; for this reason we shall retain all such diagrams here. On the other hand 
we follow Ref. || and keep only those NN- and 7r AW-irreducible connected diagrams in the 
process 7rAW — > irNN (the three-body forces) that are necessary to avoid undercounting. The 
formulation retaining all three-body forces will be given elsewhere. 

It is easy to rearrange the four- dimensional irNN equations of Ref. M into a convenient 



form similar to the one used by Afnan and Blankleider [15| in a three-dimensional formulation 



of the tiNN system. For the distinguishable nucleon case we obtain 

T d = V d + V d gfT d (1) 
where T d , V d , and Qf are 4x4 matrices given by 
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Note that we use a notation where the superscript d denotes the distinguishable nucleon case 
on those symbols that will be used later without the superscript for indistinguishable nucleons. 
Eq. ([!]) is a symbolic equation representing a Bethe-Salpeter integral equation to be solved for 
T d . To clarify its meaning we give the explicit numerical form of the equation for NN scattering 
at the end of this subsection. T d consists of transition amplitudes T$ N , T\n, Tn^, and T\ M 
(A and /i are spectator-subsystem channel labels) the last three being elements of the matrices 
Tff, ff,, and T d , respectively The physical amplitudes for NN — > NN, NN — > ird, ird^NN, 
and ird — > ird are then given by 

XffN = i X dN = ^dT 3 N ; Xff d = T N3 ifj d ; Xj d = ^ d T 3 ^ d , (3) 

respectively, where ip d is the deuteron wave function in the presence of a spectator pion. 

The elements making up the kernel V d , specified in Eq. (^), consist of the quantities V$ N , 
T d , T d , and X defined as follows: 

V$ N = V^ Ed + V^ d + F ( d G (F 1 + F 2 ) + (Ft + F 2 )G F d 

+ B d G B d + F d G F d - A (4) 

where V^ Ed is the nucleon-nucleon one pion exchange potential, vj^ d is the 7riViV-irreducible 
part of the NN potential, Fi = fidj 1 where i, j = 1,2 with i ^ j, fa being the iV, — > 7riVj vertex 
function and dj the Feynman propagator of nucleon j, F d is the simultaneously NN- and tiNN- 
irreducible connected amplitude for NN — > ttNN, B d = B + PBP where B = V^ Ed d 2 f 2 and 
P is the nucleon exchange operator, and A is a subtraction term that eliminates overcounting. 
Diagrams illustrating V^ Ed , V^ d , F 2 , B and F d are given in Fig. ||. T d is a 3 x 1 matrix 
whose A'th row element is given by 

H = E SMFi + F d - B d (5) 

where 5\i = l — 8xi- Note that here B d plays the role of a subtraction term. T d is the 1x3 matrix 
that is the time reversed version of T d (similarly for other "barred" quantities), Go is the 7riV7V 
propagator, and X is the matrix whose (A,/i)'th element is 5\^. Finally the propagator term 
Qf is a diagonal matrix consisting of the NN propagator Dq, and the 3x3 diagonal matrix 
w° whose diagonal elements are tid 2 l , t 2 di 1 , and t d d 3 l , with t\ being the two-body t matrix 
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FIG. 2. Illustration of quantities making up V^ N . (a) V^ Ed , (b) a contribution to vj^ d , (c) 
F 2 = /2^x X , (d) B, and (e) a contribution to F d . 
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FIG. 3. Terms making up the subtraction term A. (a) W ww , (b) W w n, (c) Waw, (d) X, (e) 
y, and (f) a contribution to i^- The dark circles represent the following two-body amplitudes: (a) 
full 7T7T t-matrix, (b) one-nucleon irreducible ttN t-matrix, and (c) full NN t-matrix minus the NN 
one-pion-exchange potential. 

in channel A (for A = 1 or 2, t\ is defined to be the nN t matrix with the nucleon pole term 
removed). The subtraction term A is defined with the help of Fig. [| as follows: 

A = WW + W$ N + W NN + X + Y d 

+ B d G F d + F d G B d + F 2 ^hh ND + khG^F^ _ (6) 



ND 



where F%^ is the simultaneously NN- and 7rAW-irreducible connected amplitude for NN — > 
nnNN, Gf_ is the tcttNN propagator, W d N = W nN + PW nN P and Y d = Y + PYP. The 
expression fif2Gl 7T F 2lT in Eq. (|6|) is a subtraction term for the overcounted contributions in 
(Fi+F2)GoF d in Eq. (|J). It consists of two types of contributions, one with the two intermediate 
state pions crossing, and one with them not crossing, as illustrated in Fig. [|. As our equations 
are derived with the ttNN vertex being dressed from the beginning 0, terms that contribute 
to a further dressing of the nNN vertex, like that of Fig. |](b), need to be suppressed. This 
suppression is indicated in Eq. (||) by the "no dressing" symbol |nd- Similar comments of course 
apply to the term F 27r G^ fif 2 in Eq. (Bp. 

To illustrate the numerical form of our equations we consider the case of the AW scattering 
amplitude Tfj N as given by Eq. (P and Eq. (§) 

3 

T N n = v nn + v nn d oT nn + ^2 J 7 x G w xx G T XN . (7) 

A=l 

For the purposes of the illustration we may simplify this equation by neglecting F d (but not 
B d as it is a subtraction term). In this case we obtain the following numerical integral equation 
iorT d N : 



(a) ,\! (b) 



i 



FIG. 4. Example of diagrams contributing to fif2GQ n F 2w in Eq. (^|). (a) Contribution with 
crossed pions, (b) contribution with pions not crossing. 
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FIG. 5. Graphical representation of the integral equation, Eq. (||), for the NN scattering ampli- 
tude. The top line involves the amplitude T^ N , the second line involves Tijv, the third line involves 
T2N and the last two lines involve T 3 jy. 



T nn (P1P2 1 P1P2) = V$ N (p[p 2 , P1P2) 
d A ki d 4 k 2 



d 4 h 



(2tt) 4 (2tt) 



(2tt) 4 

Mp 2 , k 2 k' 3 )d^{k 3 )t x {p' x k' 3 , k ih) 
d 4 k[ 



Vnn (p'xP'2 ,kxk 2 )Do{kik 2 ) T^ N ( fc x k 2 , Pip 2 ) 



(2tt)< 



fM, k'lk' 3 )DMKK)vt£ E (Kp 2 , k[k 2 ) 
+ f 2 ( P ' 2 , k'^DMKKW^ip'A, k[k 2 ) 



x dtik'Jhik'^, hk 3 ) I Goihhh^ihhk^p.p,) + I ^ 
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d A ki d A k 2 
(2vr) 4 (2tt) 4 



h^k'^x{K) 



t z {k' x p' 2 ,k x k 2 



d A k'[ 



V^ E (k' lP ' 2 , k'!k 2 )D (k'{k' 2 )t 3 {k'!k' 2 , hk 2 ) 



x G (hk 2 k 3 )T 3N (k 1 k 2 k 3 , p t p 2 ) + f ^ 



P2 
P2 



where to save on notation we have suppressed spin and isospin labels and used the sym- 



bol 



to indicate a term derived from the preceding contribution by exchang- 



er «-> P' 2 
px <-> p 2 

ing the momenta, spin, and isospin labels of the external nucleons. In writing Eq. (H) we 
have used the fact that numerically fi(p',pk) = f 2 (p',pk), ti(p'k',pk) = t 2 (p'k',pk), and 
T 1N (k 2 kik 3 ,P2Pi) = ^{kxhk^pxp-z). In Eq. (|) D 7tN (kik 3 ) = d^k^d^k^) is the %N propa- 
gator, and all unspecified momenta are determined by momentum conservation, e.g. in all the 
above terms k 3 = P i + p 2 — ki — k 2 . Eq. (|8|) is illustrated in Fig. ||. In reference to this figure 
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we note that despite appearances, all the subtraction terms (terms entering with a minus sign) 
end on the left in the same way, either with B = f2d 2 V^ E or with PBP = fidiV^ E . Also, 
the last two lines clearly involve the contribution t 3 — V^ E D t 3 which may appear incorrect 
as the term subtracted involves an overcounted NN interaction. Yet just such an overcounted 
contribution needs to be subtracted to stop overcounting in the overall equations. In this re- 
spect it is interesting to note that in the one pion exchange ladder approximation to t 3 , the 
combination t 3 — V^ E D t 3 reduces down simply to V^" E . 



B. Indistinguishable nucleon case 



In approaches based on second quantisation in quantum mechanics it is usual to obtain 
the scattering equations for identical particles by explicitly symmetrizing the equations of the 
distinguishable particle case. This procedure is not justified in the framework of relativistic 
quantum field theory. Nevertheless, as we have already derived the ttNN equations taking 
into account identical particle symmetry right from the beginning |J, we can formally deduce 
how the above distinguishable nucleon nNN equations need to be modified in order to get the 
indistinguishable nucleon case. With this in mind, we introduce identical nucleon transition 
amplitudes defined in terms of distinguishable nucleon transition amplitudes as 



(9) 
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T 13 A 
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= T 33 A 



where A = 1 — P is the antisymmetrizing operator. As T N1 = PT N2 P and 2^3 = PTnzP) we 
can alternatively write T N a = AT N i, T Nd = AT N3 , and T d A = AT 3 i. Thus in the transition 
to the indistinguishable particle case, the original 16 transition amplitudes for distinguishable 
particles have been reduced to 9 antisymmetric transition amplitudes. By taking residues of 
the nNN Green function for identical nucleons at two-body subsystem poles, one obtains the 
following expressions for the physical amplitudes: 

Aattv = Tnn ] XdN = IpdTdN ] ^Nd = TNdlpd ; Xdd = ^dTdd^d- (10) 

Using Eqs. (0) in Eq. ([]]) it is easy to show that one again obtains a Bethe-Salpeter equation 

T = V + Vg t T (11) 
but where now T, V, and Q t are 3x3 matrices given by 
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\Got 3 d 3 1 Gc 



f V NN A ATt HA 
TfA -Gq X P Gq 1 A 
\ T 3 A Gq 1 A 



(12) 
(13) 



where t 3 = t 3 A is the t matrix for two identical nucleons in the presence of a spectator pion. 
From Eq. (H) it follows that 
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VnnA — Vnn E + ^aw + F c GqF r + FiGqF c 

+ F^G^hh + hhGfF* - W R - W™ - W R N -X R - Y LR 

+ l -(B LR -F c )G (B LR -F c ) (14) 

where V^ E = V$$ Ed A, vj$ = V^ d A, F c = F*A, F c = F d A } and where we have used a 
superscript L to indicate that A is acting on the left, and a superscript R to indicate that A is 
acting on the right. Note that W R n = W£ w , W§ N = W^ N , X R = X L , and F£ = F 2 L n . Here we 
have also used the simple result 

(B d - F c d )G (B d - F d )A = ^(B LR - F C )G (B LR - F c ). (15) 

The other terms in the kernel V of Eq. ([12]) are 

T d = F 2 -B d + F d - T d = F 2 -B d + F d , (16) 

Tl = F l + F 2 -B d + F d - F d = F 1 + F 2 -B d + F d . (17) 

The %NN equations as given by Eqs. (|TT^ , (O) and (|T3| ) have a form where seven of the 
eight non-zero elements making up the kernel V contain the antisymmetrization operator A. 
One element (£3) in Q t also contains A. The question arises if one can remove the operators A 
from these equations by analogy with the familiar example of the one-channel Bethe-Salpeter 
equation for two identical nucleons: 

T = VA(l + ^GT). (18) 

Here the NN potential V is multiplied on the right by A, but as V has the property VP = PV 
and therefore V A = AV, it is easy to see that we can instead solve the equation 

T = V(l + GT') (19) 



where A has been removed from the equation and where the NN t matrix T is obtained by 
antisymmetrizing the solution: T = T'A. 

Unfortunately it is not possible to follow a similar procedure to completely remove the 
operator A from Eqs. (|Tl~D, (|l2j) and ([13]). The reason is that VP 7^ PV for the kernel V 
of Eq. (|I2|) (in particular J-'fP 7^ PTf). In this respect we note that the ttNN equations for 
identical particles given in Ref. are not equivalent to ours as they do not involve the operator 
A. Although A cannot be removed completely from the identical particle ttNN equations, one 
can reduce the number of places where A appears. For example, it is easy to show that Eqs. ([11]), 
( |j"2|) and (|13|) are equivalent to the equations 

T = V' + VG'tT; (20) 



V 



/ v d 

I V NN 

TfA 



~Gq P Gq a 
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Q' t =\ Gohd^Go 

Gotjd^Go 



(21) 



where 
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and all but two operators A have been removed form the integral equations. There are further 
forms of the ttNN equations with operators A appearing at various other places in the equations, 
yet there are always at least two A operators present. 



III. tvNN ELECTROMAGNETIC TRANSITION CURRENTS 
A. Gauging the ttNN equations 

In this section we shall derive expressions for the various electromagnetic transition currents 
of the nNN system. To do this we utilise the recently introduced gauging of equations method 
[0,0. As the gauging procedure is identical for the distinguishable and indistinguishable par- 
ticle cases, we restrict our attention to the tcNN system where the nucleons are treated as 
indistinguishable particles. 

As discussed in the previous section, the strong interaction ttNN equations can be written 
in a number of equivalent forms. Choosing the form given by Eqs. (0), (|12|) and fll3D , direct 
gauging of Eq. flll]) gives 

= V M + V^Q t T + Vg?T + V&T" (23) 

which can easily be solved for T M giving 

T" = (1 + Tg t )V^{l + g t T) + TQ^T. (24) 

T M is a matrix of gauged transition amplitudes T NN , T^ A , T^ d , etc. To obtain the physical 
electromagnetic transition currents of the ttNN system where photons are attached everywhere 
it is not sufficient to just gauge the physical ttNN amplitudes of Eq. (|10D . Although this 
would indeed attach photons everywhere inside the strong interaction diagrams, it would miss 
those contributions to the physical electromagnetic transition currents where the photons are 
attached to the external (initial and final state) pions and nucleons. In order to also include 
these external leg contributions it is useful to attach the corresponding propagators to the 
A-amplitudes of Eq. (0): 

Xnn = DqXmnDq XdN = d n XdND (25) 
XNd = D X Nd d n X dd = d^Xdd d n . (26) 

The physical electromagnetic transition currents are then obtained by gauging Eqs. (B^) and 
(p6|) and "chopping off" external legs: 

3nn = D 1 X I ^ }N D Q 1 j dN = d^X^D^ 1 (27) 
f Nd = D^X^ d- 1 f dd = d- l X^ d-\ (28) 

Using Eqs. (0) we obtain that 
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Jnn — D l D^T NN + T nn DqD 1 + Tfc N (29) 
3dN = 4>d D o T dN + d~ l 4> d G%T dN + (fi d D T^ N 

+ faDoTwDSDo 1 (30) 
f m = T Nd D Q <% + T Nd G& d d? + T£ d D cfi d 

+ D^D!^r Ni Do<t> d (31) 
3dd = fidDoTddDofa + d~ l (p d D^T dd DQ(j) d 

+ 4> d D T^ d D Q (f) d + faDoTddD^dd' 1 

+ faDoTuDoffi (32) 

where (fi d is the deuteron bound state vertex function defined by the relation ip d = did 2 (fi d . The 
above equations express the physical electromagnetic transition currents j£p (a,/3 = N or d) in 
terms of the half-on-shell 7riV7V transition amplitudes T a p and the gauged quantities V£p, (Gt) a p, 
and (fid- Note that (fid consists of contributions where the photon is attached everywhere inside 
the deuteron bound state, and is determined by gauging the two-nucleon bound state equation 

farina- 

That the above equations are gauge invariant is evident from the fact that we have formally 
attached the photon to all possible places in the strong interaction model. The gauge invari- 
ance of our equations also follows from a strict mathematical proof; however, as this proof is 
essentially identical to the one given for the NNN system 0, we shall not repeat it here. 



B. Alternative form of the 7rNN equations 



Although the preceding discussion solves the problem of gauging the nNN system, the 
expression obtained for the gauged transition amplitudes, Eq. (pip, may not be the most con- 
venient for numerical calculations. The disadvantage of Eq. ([24]) is that it utilises a Green 
function Q t which contains two-body t matrices, while such t matrices are already implicitly 
present in the adjoining amplitudes of T. This makes the calculation of T M unnecessarily com- 
plicated. One could eliminate this multiple t-dependence in Eq. (|4]) by making use of Eq. (|TJ); 
however, this would be a lengthy and awkward procedure. Instead we derive an alternative 
form of the ttNN equations which uses a "free" Green function which contains no two-body 
interactions and which leads to simpler expressions for the nNN electromagnetic transition 
currents. 

The ttNN equations Eqs. fllT]), ( |T2"D and (|H|), can be written in the form 



Tnn T n 
T N T 




LG Q 1 



1 



§F>o 

GntG 




(33) 



where Vnn 
T = 



V$ N A is given by Eq. flTJ) , and 
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(34) 



(35) 
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-P A 
A 



t 







(36) 



With the view of gauging Green function versions of ttNN transition amplitudes, we introduce 
the Green function matrix appropriate to Eq. 



Gnn 


Gn 


Gn 
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AD 


°) 





) 



( Gnn Gna Gn<i 
Gan Caa G&d 
\ GdN GdA Gdd 



D 
Go 



Tnn T n 
T n T 



Do 
Go 



(37) 



The inhomogeneous term is chosen so that Gnn corresponds exactly to the full Green function 
for NN scattering. Then, as shown in the Appendix, Q satisfies the equation 



AD 
LGo 



1 + 



\V m -\KLG Q K |A 



t 



where A and A are defined by Eqs. ([75|) . The essential feature of Eq. 
in terms of an effective "free" Green function matrix 



(3? 



is that it is written 



Go 



ADq 
LGo 



(39) 



which does not involve two-body interactions. For this reason Eq. (j38l) is ideal for the purposes 
of gauging. 



C. Gauging the alternative form of the wNN equations 

In terms of the elements of Q, the Green function versions of the physical amplitudes [defined 
in Eq. (pSDI are given by 

X N n = Gnn ! X d N = 4><iGdN X N d = GNd4>d ] X dd = ^dQdd^d- (40) 
After gauging, these equations give 

X NN = G NN (41) 

X&r = €G dN + 4>A (42) 

X Nd = G Nd <Pd + GNd€ (43) 

X%d = €Gdd4>d + faG^fa + 4> d G dd r d - (44) 

The ttNN electromagnetic transition currents j^p are then determined by Eqs. ([27]) and 
Eqs. (H|). 

To determine the quantities G^g in Eqs. d4l|) -(|44) we need to derive the expression for £/ M 
by gauging Eq. fl38|). Defining 



\V m -\KLG^ ±A 



V t =( * VNN if^ 1 ^ (45) 
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Eq. (|38| ) can be written as 

G = Go + G V t G. (46) 
Gauging this equation and solving for G^ gives 

Q* = (1 + 0V t )0£(l + V t G) + (47) 
To simplify this equation we cannot use Eq. ( |46|) to write 1 + VtG = Go~ X G and 1 + GVt = GGo 



-i 

since L is singular so that the inverse Gq 1 does not exist. Instead we use the fact that L = LQL, 
where 

/ i+P i \ 



which allows us to write 



where 



1 4 J (48) 



GS = GoM^Go (49) 



M* = ( i D ^ D ^ n ° n _ x ) . (50) 



QGq Gq Gq 

Using Eq. ( fj9|) in Eq. (f|7[) allows us to write a compact expression for G 

g» = g(M" + V?)Q. (51) 



Comparing this expression with the one of Eq. fl24|) , we see that both involve the gauged two- 
body t matrix t^; however, in contrast to Eq. (p4|), the above expression does not contain 
adjoining strong interaction t matrices and may therefore be preferable for numerical calcula- 
tions. 

As Dq = {did 2 )^ = d^d 2 + d\d 2 and Gq = {did 2 d^Y = d^d 2 d^ + d\d 2 d 3 + did 2 d^, we obtain 

Dq'D^Dq 1 = T^d, 1 + d?T% (52) 

Gq ^GqGq ^ — J^^d 2 ^£^3 ^ ~\~ d-y ^T^g^ ^ "I - d^ ^d 2 ^1^3 (^3) 

where 

rf = d-'d^ (54) 

is the electromagnetic vertex function of particle i. The term of Eq. (|5l|) thus corresponds 
to photon coupling in the impulse approximation. The gauged matrix Vf corresponds to the 
interaction currents and consists of the elements = (Vjvjv — ALGqA)^, A m , A m and P. It is 
important to note that the diagonal elements of matrix f are both of the form 

= tfdj 1 + ^ (^7 1 ) M = tfdj 1 - ur* (55) 

where the last equality follows from the fact that (d~j l d^ = 0. Thus the diagonal elements 
of t M involve new subtraction terms tiT^ whose origin does not lie in the subtraction terms of 
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V 



A — 



r-x- 











— A 



V 







TT 



XT 



J 



(a) (b) 
FIG. 6. Graphical representation of (a) A4^ and (b) which enter into the expression for the 
gauged Green function matrix of Eq. ( |5ll) . Constant factors and effects of operator A and matrix 
Q have been ignored in this illustration. 



the strong interaction ttNN equations, but rather in the gauging procedure itself. Analogous 
subtraction terms arise in the three-nucleon problem whose strong interaction equations have 
no subtraction terms 0,0. Similar subtraction terms will arise in the gauging of F\ and F\ 
contained in A and A respectively. The graphical representation of Ai^ and Vf is given in 
Fig- 1- 

Eq. (O) can be used to determine all the possible electromagnetic transition current of the 
ttNN system. For example, current j^ N for the electromagnetic transition NN — > NN is given 
by 

3nn = D l G tl NN D Q 1 (56) 
= £ DZ x G Na (A4 M + Vt)ap GpnD, 1 

a/3 

= J2 (AS Na + T Na D a ) (M» + V?) a(3 (A5 m + D?T m ) (57) 

a/3 

where a, (3 = N, A, d, and = D , = Dd = Go. 

An especially interesting use of Eq. ([y]) is to study the electromagnetic properties of bound 
states of the irNN system. It is certainly expected that the strong interaction irNN model 
under discussion admits a bound state corresponding to the physical deuteron. In this case a 
solution will exist to the homogeneous version of Eq. flTT|): 

$ = VQ t § (58) 

where $ is a matrix of deuteron vertex functions 



$ = $a • (59) 




Here $at is the usual deuteron vertex function describing the d — *> NN transition, while $a and 
$d are somewhat unusual in that they describe transitions to clustered nNN states: d — ► (ttN)N 
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and d — ► (NN)n respectively. Comparing Eq. ( |3?| ) and Eq. ( p8| ) it is seen that Q has a pole at 
the deuteron mass M^. 

where P is the total four-momentum of the system, and where ^ satisfies the equation 

* = £oV t ^. (61) 
Clearly \l/ is related to $ by the equation 

so that either of the equations (|58D or (^l|) can be used to determine Taking the left 
and right residues at the deuteron bound state poles of Eq. ( |5l~D we obtain the bound state 
electromagnetic current 

f = ^ {M* + Vf ) ^ (63) 

which describes the electromagnetic properties of the deuteron whose internal structure is 
described by the present nNN model. Eq. (^) provides a very rich description of the internal 
electromagnetic structure of the deuteron with all possible meson exchange currents being taken 
into account in a gauge invariant way. In view of the accuracy of this model which is based 
on meson and baryon degrees of freedom, a comparison of the deuteron electromagnetic form 
factors (easily extracted from j M ) with those extracted from experiment should prove to be 
most interesting. 



IV. CONCLUSIONS 

We have derived gauge invariant expressions for the electromagnetic transition currents 
of the nNN system where the strong interactions are described by four- dimensional integral 
equations. The feature of our approach is that the external photon is coupled everywhere in 
the strong interaction model, in this way giving a unified description of the nNN and jnNN 
systems. This unified approach to the nNN system and its electromagnetic currents has been 
made possible by the recent introduction of the gauging of equations method Jl]||. The use 
of this method has also enabled us to avoid all the overcounting problems that are inherent in 
four-dimensional descriptions of 7riV7V-like systems ||. 

The expressions we have derived can be used directly to make four-dimensional calculations 
of all the reactions induced in the nNN system by an external electromagnetic probe. However, 
in view of the practical difficulty of solving four- dimensional integral equations, it may also be 
useful to have a gauge invariant three-dimensional description of the same processes. In that 
case our four- dimensional expressions can be used to provide the starting point for a three- 
dimensional reduction. One such three-dimensional reduction scheme that preserves gauge 
invariance and is easily applied to our four- dimensional expressions was discussed in Ref. ||. 

Although we have specifically gauged the nNN system, it should be noted that the derived 
expressions apply equally well to other systems consisting of two fermions and one boson which 
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can be absorbed by the fermions. For example one could apply our expressions to the quark- 
antiquark-gluon system in order to calculate the meson spectrum including its electromagnetic 
properties. 

It is also worth noting that our derivation of the gauged irNN equations only assumed that 
the external field couples to hadrons only to first order in the field-hadron coupling constant, 
but otherwise does not depend on the nature of the external field. Thus, for example, our 
expressions can be used directly to determine the weak interaction transition currents of the 
nNN system. 
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APPENDIX 



In this Appendix we derive Eq. (j38|) . Writing Eq. fl33f) as T = V + VQ{T and Eq. (|37|) as 
Q = Qu + gTg where 

^=( A 0°l)i 9=1? I. (64) 




V 



Vj 



NN 



T LG 



-i 



Gt 




(65) 



we have that 



Q = Qu + gTg = Q u + gVg + gVQ t T g = g u + gVg + gVQ t g-\Q - Q u ) (66) 

or 

Q = Q u + gV(g- Qt.g^Qu) + gVQ t g- l Q. (67) 

Using the explicit matrix forms for Q u , g, V, and Q t given in Eq. (|64]) and Eq. (|65"D , the equation 
for Q takes the form 

r ( AD D fG \ ( \D Q V NN D fG t \ p 

G= { LG ) + [ 2 \G,T LG t ) Q - ^ 
Multiplying this equation from the left by the matrix 



where 



Z = ^ Y _^ j so that ZL = LZ = fj | j , (70) 
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we obtain 



1 -D FZ 
ZL 



AD 
LG 



D 

LGq 



\Vnn ~ \TGqZJ- 
\ZT 



t 



Q. 



(71) 



Transferring the off-diagonal term —D^TZ from the left hand side to the right and recognising 
that 



1 
LZ 



Q = Q 



(72) 



Eq. Q7ip becomes 

Q 

This equation may be written in the form 



AD 


° V 







{ o 


LGq i 







2^iVJV — ^J'GqZJ- J- Z 
\ZT t 



Q. 



(73) 



AD 
LG 



1 + 



NN 



(74) 



where 
A 



F x - \(B d -F a c 



A 



A = [A [F x - \{B d - F d )] jrf)}. (75) 



We note that the presence of the antisymmetrization operator A in the term ADq of Eq. ([74]) 
allows one to eliminate the antisymmetrization operators in Vnn and A by making the replace- 
ments 



V NN -> 2V^ and A 
in Eq. ( [74] ) and Eq. (|75|) respectively. 



f5 a 



(76) 
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